Expression Localization of the KRT32 Gene and Its Association of Genetic Variation with Wool Traits.

Changes in keratin gene expression and spatiotemporal regulation determine the compositional content and cellular localization of wool keratin, thereby affecting wool traits. Therefore, keratin gene family member 32 (KRT32) was selected for a study using RT-qPCR, immunofluorescence, and penta-primer amplification refractory mutation system (PARMS) techniques. The results showed that KRT32 mRNA was highly expressed in the skin and localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP). Sequencing results identified eight SNPs in KRT32, and association analyses revealed that the variations were significantly associated with multiple traits in wool (p < 0.05), including MFD, CF and MFC. The constructed haplotype combination H2H3 has higher CF and smaller MFD than other haplotype combination (p < 0.05). In conclusion, KRT32 can be used as a candidate gene for molecular genetic improvement of wool in Gansu Alpine Fine-wool sheep.

Molecular Characteristics of JAK2 and Its Effect on the Milk Fat and Casein Synthesis of Ovine Mammary Epithelial Cells.

In addition to its association with milk protein synthesis via the Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway, JAK2 also affects milk fat synthesis. However, to date, there have been no reports on the effect of JAK2 on ovine mammary epithelial cells (OMECs), which directly determine milk yield and milk contents. In this study, the coding sequence (CDS) region of ovine JAK2 was cloned and identified and its tissue expression and localization in ovine mammary glands, as well as its effects on the viability, proliferation, and milk fat and casein levels of OMECs, were also investigated. The CDS region of ovine JAK2, 3399 bp in length, was cloned and its authenticity was validated by analyzing its sequence similarity with JAK2 sequences from other animal species using a phylogenetic tree. JAK2 was found to be expressed in six ovine tissues, with the highest expression being in the mammary gland. Over-expressed JAK2 and three groups of JAK2 interference sequences were successfully transfected into OMECs identified by immunofluorescence staining. When compared with the negative control (NC) group, the viability of OMECs was increased by 90.1% in the pcDNA3.1-JAK2 group. The over-expression of JAK2 also increased the number and ratio of EdU-labeled positive OMECs, as well as the expression levels of three cell proliferation marker genes. These findings show that JAK2 promotes the viability and proliferation of OMECs. Meanwhile, the triglyceride content in the over-expressed JAK2 group was 2.9-fold higher than the controls and the expression levels of four milk fat synthesis marker genes were also increased. These results indicate that JAK2 promotes milk fat synthesis. Over-expressed JAK2 significantly up-regulated the expression levels of casein alpha s2 (CSN1S2), casein beta (CSN2), and casein kappa (CSN3) but down-regulated casein alpha s1 (CSN1S1) expression. In contrast, small interfered JAK2 had the opposite effect to JAK2 over-expression on the viability, proliferation, and milk fat and milk protein synthesis of OMECs. In summary, these results demonstrate that JAK2 promotes the viability, proliferation, and milk fat synthesis of OMECs in addition to regulating casein expression in these cells. This study contributes to a better comprehension of the role of JAK2 in the lactation performance of sheep.

Screening and identification of lncRNAs in preadipocyte differentiation in sheep.

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-ß) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Identification and screening of circular RNAs during adipogenic differentiation of ovine preadipocyte by RNA-seq.

Circular RNAs (circRNAs) are a class of non-coding RNAs that play important roles in preadipocyte differentiation and adipogenesis. However, little is known about genome-wide identification, expression profile, and function of circRNAs in sheep. To investigate the role of circRNAs during ovine adipogenic differentiation, the subcutaneous adipose tissue of Tibetan rams was collected in June 2022. Subsequently, the preadipocytes were immediately isolated from collected adipose tissue and then induced to begin differentiation. The adipocytes samples cultured on days 0, 2, and 8 of preadipocytes differentiation were used to perform RNA sequencing (RNA-seq) analysis to construct the expression profiles of circRNAs. Subsequently, the function of differentially expressed circRNAs was investigated by performing the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of their parent genes. Finally, a circRNAs-miRNAs-mRNAs network involved in adipogenic differentiation was been analyzed. As a result, a total of 6,449 candidate circRNAs were identified in ovine preadipocytes. Of these circRNAs identified, 63 candidate circRNAs were differentially expressed among the three differentiation stages and their parent genes were mainly enriched in acetyl-CoA metabolic process, positive regulation of lipid biosynthetic process, positive regulation of steroid biosynthetic process, and focal adhesion pathway (P < 0.05). Based on a circRNAs-miRNAs-mRNAs regulatory network constructed, circ_004977, circ_006132 and circ_003788 were found to function as competing endogenous RNAs (ceRNAs) to regulate ovine preadipocyte differentiation and lipid metabolism. The results provide an improved understanding of functions and molecular mechanisms of circRNAs underlying ovine adipogenesis in sheep.

The moderate fat deposition contributes to improve mutton quality, which is associated with the differentiation of preadipocytes. To investigate roles of circular RNAs (circRNAs) in preadipocyte differentiation, we identified circRNAs on days 0, 2, and 8 of preadipocytes differentiation and compared the expression profile of circRNAs at different adipogenic differentiation stages. A total of 6,449 candidate circRNAs were identified, among which 63 candidate circRNAs were differentially expressed among the three differentiation stages. The parent genes of differentially expressed circRNAs were enriched in several biological process and pathways related to lipid metabolism and synthesis. In addition, several circRNAs may regulate ovine preadipocyte differentiation by interacting with microRNAs (miRNAs). The results reveal the potential roles of circRNAs in adipogenic differentiation of sheep.

Spatiotemporal Expression and Haplotypes Identification of KRT84 Gene and Their Association with Wool Traits in Gansu Alpine Fine-Wool Sheep.

Keratin (K) is a major protein component of hair and is involved in hair growth and development. In this study, we analysed the expression, localization, and polymorphism of the K84 gene (KRT84) in Gansu Alpine Fine-wool sheep using immunofluorescence, RT-qPCR, and PARMS (penta-primer amplification refractory mutation system). Haplotypes of KRT84 were also constructed and their relationship with wool traits analysed. It was revealed that KRT84 was highly expressed in hair follicles, including the inner root sheath, outer root sheath, and hair medulla and at all six lamb ages investigated from 1 to 270 days of age. Three SNPs were detected in KRT84 exon 1, and they formed three haplotypes (named H1, H2, and H3) and six genotypes. Analyses revealed an association between haplotype combinations (diplotypes) and the mean fibre curvature, mean staple length, mean staple strength, mean fibre diameter, the coefficient of variation of fibre diameter, and comfort factor for these sheep. These results suggest that KRT84 is of importance in determining several key traits in Gansu Alpine Fine-wool sheep and that the gene could possibly be used as a genetic marker for wool trait selection in these sheep.

Synergistic Responses of Tibetan Sheep Rumen Microbiota, Metabolites, and the Host to the Plateau Environment.

Plateau adaptation in animals involves genetic mechanisms as well as coevolutionary mechanisms of the microbiota and metabolome of the animal. Therefore, the characteristics of the rumen microbiome and metabolome, transcriptome, and serum metabolome of Tibetan sheep at different altitudes (4500 m, 3500 m, and 2500 m) were analyzed. The results showed that the rumen differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and carbohydrate metabolism, and there was a significant correlation with microbiota. The differentially expressed genes and metabolites at middle and high altitudes were coenriched in asthma, arachidonic acid metabolism, and butanoate and propanoate metabolism. In addition, the serum differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and metabolism of xenobiotics by cytochrome P450, and they were also related to microbiota. Further analysis revealed that rumen metabolites accounted for 7.65% of serum metabolites. These common metabolites were mainly enriched in metabolic pathways and were significantly correlated with host genes (p < 0.05). This study found that microbiota, metabolites, and epithelial genes were coenriched in pathways related to lipid metabolism, energy metabolism, and immune metabolism, which may be involved in the regulation of Tibetan sheep adaptation to plateau environmental changes.

Functional differentiation of the ovine preadipocytes -insights from gene expression profiling.

The preadipocytes differentiation is a vital process of lipogenesis; exploring the molecular mechanisms of lipogenesis contributes to improve the meat quality and final commercial income. Lipogenesis has been widely reported in other livestock, but little is known about the gene expression profiles at different stages during preadipocytes differentiation in sheep. In this study, ovine preadipocytes were cultured in vitro and then induced to begin differentiation. Then, the gene expression profiles of preadipocytes collected on day 0 (D0), day 2 (D2), and day 8 (D8) of differentiation were analyzed by RNA-seq technology. According to the findings, 2254 differentially expressed genes (DEGs) were found in D2 vs D0; 1817 DEGs and 1902 DEGs were found in D8 vs D0 and D8 vs D2, respectively. The DEGs were found to be enriched in several biological processes, including focal adhesion, ECM-receptor interaction, PI3K-Akt signaling pathway, steroid biosynthesis, and MAPK signaling pathway, according to Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The regulatory network of the DEGs related to ovine preadipocytes differentiation was systematically constructed, which showed that hub genes might modulate ovine preadipocytes differentiation. In summary, preadipocyte differentiation is regulated by several key genes, including ACACB, CXCL6, SREBF1, INSIG1, APOE, GJA1, CDH11, SYNE1, PCSK1, S100A4, FN1, PLIN2, CXCL6, FN1, PTX3, and FABP3. This study provides a deeper knowledge of the roles of genes in sheep lipogenesis by revealing global gene expression profiles during preadipocyte differentiation.

Identification and Characterization of Circular RNAs (circRNAs) Using RNA-Seq in Two Breeds of Cashmere Goats.

Circular RNA (circRNA) is a type of non-coding RNA generated from back-splicing the reactions of linear RNA. It plays an important role in various cellular and biological processes. However, there are few studies about the regulatory effect of circRNAs on cashmere fiber traits in cashmere goats. In this study, the expression profiles of circRNAs in skin tissue were compared between Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats, with a significant difference in cashmere fiber yield, cashmere fiber diameter, and cashmere fiber color, using RNA-seq. A total of 11,613 circRNAs were expressed in the caprine skin tissue, and their type, chromosomal distribution, and length distribution were characterized. A total of 115 up-regulated circRNAs and 146 down-regulated circRNAs in LC goats were screened compared to ZB goats. The authenticity of 10 differentially expressed circRNAs was validated by detecting their expression levels and the head-to-tail splice junction using RT-PCR and DNA sequencing, respectively. The parent genes of differentially expressed circRNA were mainly enriched in some Gene Ontology (GO) terms and pathways related to cashmere fiber traits, such as the canonical Wnt signaling pathway, which is involved in the regulation of cell promotion, stem cell proliferation, Wnt signaling pathway regulation, epithelial morphogenesis, MAPK signaling pathway, and cell adhesion molecules pathway. Eight differentially expressed circRNAs were further selected to construct a circRNA-miRNA network, and some miRNAs that were previously reported as related to fiber traits were found in the network. This study provides a deep understanding of the roles of circRNAs in the regulation of cashmere fiber traits in cashmere goats and the involvement of differential splicing in phenotypic expression according to breed and region.

RNA-Seq Reveals the Roles of Long Non-Coding RNAs (lncRNAs) in Cashmere Fiber Production Performance of Cashmere Goats in China.

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA being >200 nucleotides in length, and they are found to participate in hair follicle growth and development and wool fiber traits regulation. However, there are limited studies reporting the role of lncRNAs in cashmere fiber production in cashmere goats. In this study, Liaoning cashmere (LC) goats (n = 6) and Ziwuling black (ZB) goats (n = 6) with remarkable divergences in cashmere yield, cashmere fiber diameter, and cashmere color were selected for the construction of expression profiles of lncRNAs in skin tissue using RNA sequencing (RNA-seq). According to our previous report about the expression profiles of mRNAs originated from the same skin tissue as those used in the study, the cis and trans target genes of differentially expressed lncRNAs between the two caprine breeds were screened, resulting in a lncRNA-mRNA network. A total of 129 lncRNAs were differentially expressed in caprine skin tissue samples between LC goats and ZB goats. The presence of 2 cis target genes and 48 trans target genes for the differentially expressed lncRNAs resulted in 2 lncRNA-cis target gene pairs and 93 lncRNA-trans target gene pairs. The target genes concentrated on signaling pathways that were related to fiber follicle development, cashmere fiber diameter, and cashmere fiber color, including PPAR signaling pathway, metabolic pathways, fatty acid metabolism, fatty acid biosynthesis, tyrosine metabolism, and melanogenesis. A lncRNA-mRNA network revealed 22 lncRNA-trans target gene pairs for seven differentially expressed lncRNAs selected, of which 13 trans target genes contributed to regulation of cashmere fiber diameter, while nine trans target genes were responsible for cashmere fiber color. This study brings a clear explanation about the influences of lncRNAs over cashmere fiber traits in cashmere goats.

The Identification and Characteristics of miRNAs Related to Cashmere Fiber Traits in Skin Tissue of Cashmere Goats.

microRNAs (miRNAs) are involved in the regulation of biological phenomena by down-regulating the expression of mRNAs. In this study, Liaoning cashmere (LC) goats (n = 6) and Ziwuling black (ZB) goats (n = 6) with different cashmere fiber production performances were selected. We supposed that miRNAs are responsible for the cashmere fiber trait differences. To test the hypothesis, the expression profiles of miRNAs from the skin tissue of the two caprine breeds were compared using small RNA sequencing (RNA-seq). A total of 1293 miRNAs were expressed in the caprine skin samples, including 399 known caprine miRNAs, 691 known species-conserved miRNAs, and 203 novel miRNAs. Compared with ZB goats, 112 up-regulated miRNAs, and 32 down-regulated miRNAs were found in LC goats. The target genes of the differentially expressed miRNAs were remarkably concentrated on some terms and pathways associated with cashmere fiber performance, including binding, cell, cellular protein modification process, and Wnt, Notch, and MAPK signaling pathways. The miRNA-mRNA interaction network found that 14 miRNAs selected may contribute to cashmere fiber traits regulation by targeting functional genes associated with hair follicle activities. The results have reinforced others leading to a solid foundation for further investigation of the influences of individual miRNAs on cashmere fiber traits in cashmere goats.

Integrated transcriptome analysis reveals roles of long non-coding RNAs (lncRNAs) in caprine skeletal muscle mass and meat quality.

Long non-coding RNAs (lncRNAs) play important roles in the growth and development of skeletal muscle. However, there is limited information on goats. In this study, expression profiles of lncRNAs in Longissimus dorsi muscle from Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with divergent meat yield and meat quality were compared using RNA-sequencing. Based on our previous microRNA (miRNA) and mRNA profiles obtained from the same tissues, the target genes and binding miRNAs of differentially expressed lncRNAs were obtained. Subsequently, lncRNA-mRNA interaction networks and a ceRNA network of lncRNA-miRNA-mRNA were constructed. A total of 136 differentially expressed lncRNAs were identified between the two breeds. Fifteen cis target genes and 143 trans target genes were found for differentially expressed lncRNAs, and they were enriched in muscle contraction, muscle system process, muscle cell differentiation, and p53 signaling pathway. A total of 69 lncRNA-trans target gene pairs were constructed, with close relationship with muscle development, intramuscular fat deposition, and meat tenderness. A total of 16 lncRNA-miRNA-mRNA ceRNA pairs were identified, of which some reportedly associated with skeletal muscle development and fat deposition were found. The study will provide an improved understanding of the roles of lncRNAs in caprine meat yield and meat quality.

Rumen Epithelial Development- and Metabolism-Related Genes Regulate Their Micromorphology and VFAs Mediating Plateau Adaptability at Different Ages in Tibetan Sheep.

The rumen is an important hallmark organ of ruminants and plays an important role in the metabolism and immune barrier of Tibetan sheep on the Plateau. However, there are few studies on rumen development and metabolism regulation in Tibetan sheep at different ages. Here, we comprehensively analyzed the immune function, fermentation function, rumen epithelial micromorphology and transcriptome profile of Tibetan sheep at different ages. The results showed that the concentration of IgG decreased and the concentration of IgM increased with age (p < 0.05), and the highest concentration of IgA was observed at 1.5 and 3.5 years of age. In terms of rumen fermentation characteristics, VFAs of 4-month-old lambs were the highest, followed by VFAs and NH3-N of Tibetan sheep at 3.5 years of age. Hematoxylin-eosin staining and transmission electron microscopy section examination of rumen epithelial tissue showed that the rumen papilla width increased with age (p < 0.001), the thickness of the stratum corneum decreased, the cells in the stratum corneum showed accelerated migration and the thickness of the rumen muscle layer increased (p < 0.001). Desmosomal junctions between the layers of rumen epithelium increased at 1.5 and 3.5 years old, forming a compact barrier structure, and the basal layer had more mitochondria involved in the regulation of energy metabolism. RNA-seq analysis revealed that a total of 1006 differentially expressed genes (DEGs) were identified at four ages. The DEGs of Tibetan sheep aged 4 months and 6 years were mainly enriched in the oxidation−reduction process and ISG15-protein conjugation pathway. The 1.5 and 3.5-year-olds were mainly enriched in skeletal muscle thin filament assembly, mesenchyme migration and the tight junction pathway. WGCNA showed that DEGs related to rumen microbiota metabolite VFAs and epithelial morphology were enriched in "Metabolism of xenobiotics by cytochrome P450, PPAR signaling pathway, Butanoate metabolism pathways" and participated in the regulation of rumen epithelial immune and fermentation metabolism functions of Tibetan sheep at different ages. This study systematically revealed the regulatory mechanism of rumen epithelial development and metabolism in the plateau adaptation of Tibetan sheep, providing a new approach for the study of plateau adaptation.

MicroRNA-200c Affects Milk Fat Synthesis by Targeting PANK3 in Ovine Mammary Epithelial Cells.

Milk fat is the foremost nutrient of milk and a vital indicator in evaluating milk quality. Accumulating evidence suggests that microRNAs (miRNAs) are involved in the synthesis of milk fat. The miR-200c is closely related to lipid metabolism, but little is known about its effect on the synthesis of milk fat in MECs of ewes. Herein, the effect of miR-200c on the proliferation of ovine mammary epithelial cells (MECs) and its target relationship with a predicted target gene were investigated. The regulatory effects of miR-200c on the expression of the target genes and the content of triglycerides in ovine MECs were further analyzed. The results revealed that the expression level of miR-200c was differentially expressed in both eight tissues selected during lactation and in mammary gland tissues at different physiological periods. Overexpression of miR-200c inhibited the viability and proliferation of ovine MECs, while inhibition of miR-200c increased cell viability and promoted the proliferation of ovine MECs. Target gene prediction results indicated that miR-200c would bind the 3'UTR region of pantothenate kinase 3 (PANK3). Overexpression of miR-200c reduced the luciferase activity of PANK3, while inhibition of miR-200c increased its luciferase activity. These findings illustrated that miR-200c could directly interact with the target site of the PANK3. It was further found that overexpression of miR-200c reduced the expression levels of PANK3 and, thus, accelerated the synthesis of triglycerides. In contrary, the inhibitor of miR-200c promoted the expression of PANK3 that, thus, inhibited the synthesis of triglycerides in ovine MECs. Together, these findings revealed that miR-200c promotes the triglycerides synthesis in ovine MECs via increasing the lipid synthesis related genes expression by targeting PANK3.

MicroRNA-381 Regulates Proliferation and Differentiation of Caprine Skeletal Muscle Satellite Cells by Targeting PTEN and JAG2.

In our previous study, microRNA (miR)-381 was found to be the most down-regulated miRNA in skeletal muscle of Liaoning cashmere goats with higher skeletal muscle mass, but the molecular mechanism involved remains unclear. In this study, primary caprine skeletal muscle satellite cells (SMSCs) were isolated and identified. We investigated the effect of miR-381 on the viability, proliferation and differentiation of caprine SMSCs, and the target relationships of miR-381 with jagged canonical Notch ligand 2 (JAG2) and phosphatase and tensin homolog (PTEN). Cells isolated were positive for SMSC-specific marker protein Pax7. This suggests that purified SMSCs were obtained. The expression level of miR-381 achieved a peak value on day 4 after SMSC differentiation, and miR-381 also significantly increased the expression levels of myogenic differentiation marker genes: myosin heavy chain (MyHC), myogenin (MyoG) and myocyte enhancer factor 2C (MEF2C) in differentiated SMSCs, the area of MyHC-positive myotubes and the myogenic index. These findings suggest that miR-381 promoted myogenic differentiation of caprine SMSCs. The CCK8 assay and EDU staining analysis showed that miR-381 mimic both inhibited the viability of SMSCs and decreased the percentage of EDU-labeled positive SMSCs. In contrast, miR-381 inhibitor had the opposite effect with miR-381 mimic. A dual luciferase reporter assay verified that miR-381 can target JAG2 and PTEN by binding to the 3'-untranslated regions (3'-UTR) of the genes. The transfection of miR-381 mimic into caprine SMSCs resulted in decreases in expression levels of JAG2 and PTEN, while miR-381 inhibitor increased the two target genes in expression. This is the first study to reveal the biological mechanisms by which miR-381 regulates caprine SMSC activities.

Characteristics and Expression of circ_003628 and Its Promoted Effect on Proliferation and Differentiation of Skeletal Muscle Satellite Cells in Goats.

In our previous a study, circ_003628 was one of the most highly expressed circular RNAs (circRNAs) in the Longissimus dorsi muscle of goats found by RNA-seq, suggesting that the circRNA may be important for caprine muscle growth and development. However, there have been no reports describing the molecular mechanisms by which circ_003628 regulates the activities of goat skeletal muscle satellite cells (SMSCs). In this study, reverse transcriptase-PCR (RT-PCR) and DNA sequencing were used to validate the authenticity of circ_003628, and its characteristics, expression profile and effect on goat SMSCs were also studied using real-time quantitative-PCR (RT-qPCR), EdU, CCK-8 and immunofluorescence assays. Circ_003628 is partially originated from 13 exons, 12 introns and 3'-untranslated regions (UTR) of caprine Myosin Heavy Chain 1 (MYH1), and 25 exons and 5' UTR of Myosin Heavy Chain 4 (MYH4), as well as intergenic sequences between the two genes. A total of 77.07% of circ_003628 were located in the nuclei of goat SMSCs, while 22.93% were expressed in the cytoplasm. The circRNAs were only expressed in triceps brachii, quadriceps femoris and longissimus dorsi muscle tissues in nine caprine tissues investigated, with the highest expression level in longissimus dorsi muscle. The expression level of circ_003628 gradually increased during differentiation periods of goat SMSCs and reached the maximum on day 6 after differentiation. The small interfering RNA of circ_003628 (named si-circ_003628) inhibited the viability and proliferation of goat SMSCs, and also decreased the expression of four cell proliferation marker genes: paired box 7 (Pax7), cyclin-dependent kinase 2 (CDK2), CDK4 and CyclinD1 in goat SMSCs. Transfection of si-circ_003628 significantly decreased the area of MyHC-labeled myotubes of goat SMSCs, as well as the expression levels of three differentiation marker genes: myosin heavy chain (MyHC), myogenin (MyoG), and myocyte enhancer factor 2C (MEF2C). These results suggest that circ_003628 promotes the viability, proliferation, and differentiation of goat SMSCs, and they also provide an improved understanding of the roles of circ_003628 in skeletal muscle growth and development in goats.

MicroRNA-199a-3p regulates proliferation and milk fat synthesis of ovine mammary epithelial cells by targeting VLDLR.

In our previous study, microRNA (miR)-199a-3p was found to be the most upregulated miRNA in mammary gland tissue during the non-lactation period compared with the peak-lactation period. However, there have been no reports describing the function of miR-199a-3p in ovine mammary epithelial cells (OMECs) and the biological mechanisms by which the miRNA affects cell proliferation and milk fat synthesis in sheep. In this study, the effect of miR-199a-3p on viability, proliferation, and milk fat synthesis of OMECs was investigated, and the target relationship of the miRNA with very low-density lipoprotein receptor (VLDLR) was also verified. Transfection with a miR-199a-3p mimic increased the viability of OMECs and the number of Edu-labeled positive OMECs. In contrast, a miR-199-3p inhibitor had the opposite effect with the miR-199a-3p mimic. The expression levels of three marker genes were also regulated by both the miR-199a-3p mimic and miR-199-3p inhibitor in OMECs. Together, these results suggest that miR-199a-3p promotes the viability and proliferation of OMECs. A dual luciferase assay confirmed that miR-199a-3p can target VLDLR by binding to the 3'-untranslated regions (3'UTR) of the gene. Further studies found a negative correlation in the expression of miR-199a-3p with VLDLR. The miR-199a-3p mimic decreased the content of triglycerides, as well as the expression levels of six milk fat synthesis marker genes in OMECs, namely, lipoprotein lipase gene (LPL), acetyl-CoA carboxylase alpha gene (ACACA), fatty acid binding protein 3 gene (FABP3), CD36, stearoyl-CoA desaturase gene (SCD), and fatty acid synthase gene (FASN). The inhibition of miR-199a-3p increased the level of triglycerides and the expression of LPL, ACACA, FABP3, SCD, and FASN in OMECs. These findings suggest that miR-199a-3p inhibited milk fat synthesis of OMECs. This is the first study to reveal the molecular mechanisms by which miR-199a-3p regulates the proliferation and milk fat synthesis of OMECs in sheep.

Tissue-Specific Expression of Circ_015343 and Its Inhibitory Effect on Mammary Epithelial Cells in Sheep.

Circular RNAs (circRNAs) are a kind of non-coding RNA that have an important molecular function in mammary gland development and lactation of mammals. In our previous study, circ_015343 was found to be highly expressed in the ovine mammary gland tissue at the peak-lactation period by using RNA sequencing (RNA-seq). In the present study, the authenticity of circ_015343 was confirmed by using reverse transcriptase-polymerase chain reaction (RT-PCR) analysis and Sanger sequencing. The circ_015343 was derived from the complete 10 exons of aminoadipic semialdehyde synthase (AASS), ranging from exon 2 to exon 11 and mainly located in cytoplasm of ovine mammary epithelial cells. The circRNA was found to be expressed in eight ovine tissues, with the highest expression level in the mammary gland and the least expression in Longissimus dorsi muscle. The circ_015343 had a lower level of expression in a sheep breed with higher milk yield and milk fat content. The disturbed circ_015343 increased the viability and proliferation of the ovine mammary epithelial cells. The inhibition of circ_015343 also increased the expression levels of three milk fat synthesis marker genes: acetyl-coenzyme A carboxylase alpha (ACACA), fatty acid-binding protein 4 (FABP4), and sterol regulatory element-binding protein 1 (SREBP1), as well as three proliferation-related genes: cyclin dependent kinase 2 (CDK2), cyclin dependent kinase 4 (CDK4) and proliferating cell nuclear antigen (PCNA), but decreased the expression level of its parent gene AASS. A circRNA-miRNA-mRNA interaction network showed that circ_015343 would bind some microRNAs (miRNAs) to regulate the expression of functional genes related to the development of mammary gland and lactation. This study contributes to a better understanding of the roles of circ_015343 in the mammary gland of sheep.

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds.

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

Identification and characterization of circular RNAs in Longissimus dorsi muscle tissue from two goat breeds using RNA-Seq.

Circular RNAs (circRNAs) are a class of non-coding RNA that play crucial roles in the growth and development of skeletal muscle. However, little is known about the role of circRNAs in caprine skeletal muscle. In this study, the size of muscle fiber and the expression profiles of circRNAs were compared in Longissimus dorsi muscle of Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotypic differences in meat production performance, using hematoxylin and eosin staining and RNA-Seq, respectively. The size of muscle fiber in LC goats was larger than those in ZB goats (P < 0.05). A total of 10,875 circRNAs were identified and 214 of these were differentially expressed between the two caprine breeds. The parent genes of differentially expressed circRNAs were mainly enriched in connective tissue development, Rap1, cGMP-PKG, cAMP and Ras signaling pathway. In conclusion, circRNAs may play important roles in skeletal mass, meat production performance and meat quality traits in goats. The results provide an improved understanding of the functions of circRNAs in skeletal muscle development of goats.

Variation in caprine KRTAP1-3 and its association with cashmere fibre diameter.

Keratin-associated proteins (KAPs) are components of cashmere fibres. The gene encoding the KAP1-3 protein (KRTAP1-3) has been described in goats, but little is known about sequence variation in this gene and if it affects cashmere fibre traits. In this study, we used a polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) technique to screen for nucleotide sequence variation in caprine KRTAP1-3 in 327 Longdong cashmere goats, then analysed association between the genetic variation that was revealed and some cashmere fibre traits. Six PCR-SSCP patterns representing six different variant sequences of KRTAP1-3 (named A to F) were revealed. Among these variant sequences, seven single nucleotide polymorphisms (SNPs) were detected, with two of them being non-synonymous. Goats with genotype AC had higher mean fibre diameter (MFD) than those with genotype AB (P < 0.001), while goats with genotype AB had higher MFD than those with AA (P < 0.001). The presence of C (P < 0.001) and B (P = 0.006) in a genotype was associated with increased MFD, and together this suggests that variation in caprine KRTAP1-3 affects the key fibre trait of MFD.